The superintegrable chiral Potts quantum chain and 
generalized Chebyshev polynomials 



G. von Gehlen" and Shi-shyr Roan 6 

a Physikalisches Institut der Universitat Bonn 
Nussallee 12, 53115 Bonn, Germany 
e-mail: gehlen@th.physik. uni-bonn. de 

Institute of Mathematics, Academia Sinica 
q . Taipei, Taiwan 

' e-mail: maroan@ccvax. sinica. edu.tw 

(N : 

u ■ 

^: 

Abstract 

| Finite-dimensional representations of Onsager's algebra are characterized by the zeros 

of truncation polynomials. The Z/\r-chiral Potts quantum chain hamiltonians (of which the 
y—i , Ising chain hamiltonian is the N = 2 case) are the main known interesting representations 

^ ■ of Onsager's algebra and the corresponding polynomials have been found by Baxter and 

Albertini, McCoy and Perk in 1987-89 considering the Yang-Baxter-integrable 2-dimensional 
chiral Potts model. We study the mathematical nature of these polynomials. We find that 
\ for N > 3 and fixed charge Q these don't form classical orthogonal sets because their pure 

recursion relations have at least N + 1-terms. However, several basic properties are very 
t— I \ similar to those required for orthogonal polynomials. The N + 1-term recursions are of the 

simplest type: like for the Chebyshev polynomials the coefficients are independent of the 
| degree. We find a remarkable partial orthogonality, for N — 3, 5 with respect to Jacobi-, 

and for N — 4, 6 with respect to Chebyshev weight functions. The separation properties of 
the zeros known from orthogonal polynomials are violated only by the extreme zero at one 
end of the interval. 

> : 

1 Introduction 

Onsager's algebra (OA) has played a crucial role in the early history of exactly solvable statistical 
systems: it served as the tool for the first solution of the two-dimensional Ising model EJ . How- 
ever, since soon after many other methods for solving the Ising model have been invented, for 
many years Onsager's algebra received little attention. It reappeared when the superintegrable 
Zjv-symmetrical chiral Potts quantum chain (SCPC) was introduced as a natural generalization 
of the transverse Ising quantum chain || §], and was shown to be integrable 0] because it forms 
representations of the OA (or the closely related Dolan-Grady algebra g, [|). The actual solu- 
tion for the lowest energy eigenvalues of the SCPC by Baxter m and Albertini, McCoy and Perk 
|| did not use the representation theory of Onsager's algebra but rather inversion relations, or, 
more generally, functional relations for the transfer matrix of the 2-dimensional Yang-Baxter 
integrable chiral Potts model j|, In 1990 B.Davies constructed finite-dimensional repre- 
sentations of Onsager's algebra. In his approach zeros of truncation polynomials play a crucial 
role and the representations are determined by a set of zeros and by s/(2, C) representations. 
While still it is not known how to obtain the Onsager algebra truncation polynomials directly 



1 



from the SCPC hamiltonian, equivalent polynomials have been obtained from the 2-dimensional 
superintegrable Potts model functional relations in |7], |8|. 

Given the central role of the truncation polynomials for models solvable due to Onsager's 
algebra, it seems desirable to improve our understanding of their mathematical properties. In 
the Ising case, after a simple mapping, these are Chebyshev polynomials (a fact implicit e.g. 
in p2]]). For N > 3 we find many properties reminiscent of, but not in full agreement with 
those of orthogonal polynomials. Looking into recursion relations, differential equations, zero 
distributions and eventual weight functions for orthogonality, we find simple structures and in 
particular, a remarkable partial orthogonality with Jacobi-weight functions. 

In the next Section we review Onsager's algebra, the appearance of the polynomials and 
Baxter's results for the Zjv-polynomials. In Sec. 3 we discuss three different versions of the 
polynomials with different locations of the zeros and check how these look in the Ising case. 
Sec. 4 gives some properties of the ^-polynomials with the zeros on the negative real axis. Sec. 5 
gives our main results, which are for the polynomials with zeros in the interval — 1 < c < +1: 
recursion relations, differential equations, zeros separation and, finally, partial orthogonality. 
Sec. 6 contains our conclusions. 

2 Finite dimensional representations of Onsager's algebra 
2.1 Onsager's algebra 

Onsager's algebra A is formed [|lj] from elements A m , G[, m € Z, I € N, I > m, satisfying 

[A h A m ] =4G,_ m ; [G h A m ] = 2 A m+l - 2 A m ^; [G h G m ]=0. (1) 

Eqs.(Q) imply an infinite set of constraints: 

[^4m+i ) A m ] = [A m , A rn —i] (2) 

and the existence of the infinite set of commuting operators Q m : 

Q m = \{A m + A^ rn + k{A m+l +A_ m+1 )); [Q h Q m ] = 0. (3) 

where A; is a parameter which usually is taken to be real. Aq and A\ generate A if they satisfy 
the Dolan-Grady-relations 

[A , [Aq, [Aq, A 1 }}} = 16 [Aq, A 1 ] ; [A u [A x , [Ai,A ] } } = 16 [A u Aq] . (4) 

Finite dimensional representations of A are obtained Jl3| , [14|, [TI|] if the A m satisfy a pure finite 
recurrence or difference equation: 

n 

a k A k _i =0 (5) 

k=—n 

(implying the same equation for the Gi). For solving this relation B.Davies 
polynomial 

T{z) = «fc zk+n - (6) 

k=—n 

From A it follows that a k is either even or odd in k: at = ±a_fc, so that 

T{z) =±z 2n+1 F{l/z) (7) 



13] introduced the 
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and the zeros of J~(z) come in reciprocal pairs zj, z- 1 (j = 1, . . . , n). The A m and G m can be 
now be expressed in terms of a set of operators E^ 1 , Hf. 

n n 

A m = 2Y, (zT El + zj m Ej) ; G m = £ (zf - zj m ) Hj (8) 

which from A obey sl(2, C)-commutation rules: 

[E+,E k ] = 5 jk H k ; [H jt Eft = ±28 jk Ef. (9) 

So A is isomorphic to a subalgebra of the loop algebra of a direct sum of sl(2, C) algebras. 

We shall be interested in the eigenvalues of hermitian hamiltonians TL ~ Aq + kA\ (Tt ~ Qq 
of eq.(^[)), with parameter k and Aq and A\ satisfying (|j). Write E^ = J x j ± iJyj, then in a 

representation i?(n, s) characterized by z\,...,z n and by a spin-s representation of all the 



Jj, we have 



n 

(A Q + kA{) zM = 2j2{V + k(z j + zri))jV+i(z j -z- 1 )J$} 



i=i 



= 4 ^ ^/l + 2k Cj + k 2 J' 
where J' x j is a 5C/ (2)-operator and 



i=i 



COS0J = 3(«,- + «7 1 ) (10) 



with 9j real for hermitian "H. 



2.2 The superintegrable chiral Potts quantum chain 

The main example of hamiltonians 7i satisfying the above conditions (for generalisations see 
|]| Q) are the ZTv-SCPC-hamiltonians §, g) given by 

L N-l 

« w = - E E r4n + tz i z f«) <"> 

i=i i=i 

where and Xj are Z^r-spin operators acting in the vector spaces at the sites j (j = 
1,...,L): 

Zi Xj = Xj Zi ; Zf = Xf = 1; w = e 2 " /iV . (12) 

The normalization agrees with (g|) if we write = -^NiAo + kAx). 

In order to find the representations Z(n,s) corresponding to the various sectors of (11), no 
direct way is known. Albertini et al. § have recursively found the low-k finite-size ground 
state energy using only eq.(|l~3|) below and the order ^-approximation of (|i~l|), but this has not 
lead to a closed formula for the polynomials. However, these polynomials have been derived by 
solving the two-dimensional chiral Potts model |], 10 1, which for the "superintegrable" choice 
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of parameters^ contains the hamiltonian ( |TT|) as a transfer matrix derivative. First, using an 
inversion relation, Baxter |?J found one low-lying sector which for k — > contains the ground 
state. Then, for the Z3-case, using functional relations for the transfer matrix, Albertini, McCoy 
and Perk ||] obtained also the excited level sectors which involve solutions of Bethe-equations. 
Finally, Baxter [llj , exploiting functional relations derived from the relation of the 2-dimensional 
integrable chiral Potts model to the six- vertex model [22], solved the general ZAr-superintegrable 
case. 0. For all sectors of (11) the spin representation turns out to be s = \ or the trivial one. 
Accordingly, all eigenvalues of the hamiltonians (11) have the form 



N 



= I a + b k + 4 ^ mj Jl + 2 k cos B j + k 2 



(13) 



a and b are integers originating from the trivial representation. 

In this talk, we will consider only the low-lying (no Bethe-excitations) sector polynomials of 
Baxter |j~l| |. Baxter's polynomials are written in the variable t or s = t which is related to 
z or c = cos 9 by 



z + z 



l + s 
1 - s 



(14) 



The polynomials are (for details we refer to Baxter's paper [11], for simplicity we put there 
P b = 0): 



Q [S) ~ N ^ [l-uHj 

3=0 V 7 



(<Jty 



-Pa. 



-P a = Q + r + L mod N. 



(15) 



Here Q denotes the Z^y-charge sector (Q = 0, 1, . . . , N — 1), and r (r = 0, 1, . . . ,N— 1) labels the 
boundary condition to ( |Tl|) defined by Zl+i = oj r Z\. Since in (|l5| ) Q and r appear only in the 
sum Q + r, we can restrict us to consider the periodic case r = 0. The degree of the polynomials 
Pq (s) in the variable s is 

~(N-1)L-Q~ 



b L ,Q 



N 



(16) 



where [x] denotes the integer part of x. Inversion of the variable s leads to simple relations: 



P { q\s) = s°^P J 



N-(Q+L) 



is' 1 ), 



(17) 



where the charge index N — (Q + L) is understood mod N. These relations correspond to the 
self-reverse property (|7|) of the interpolation polynomials. 



Considering sequences of these polynomials for fixed Q and L E N, from (16) we see that 
the dimensions b^n do not always increase by one when increasing L by one: at every Nth 
step the dimension stays the same: the dimensions of the Pq for L + Q mod N = and 
L + Q mod N = 1 conincide, see e.g. Table 1. 

For N = 3 the definition (|i~5|), written explicitly, gives 



P { qIM 



-Pa 



{(*2 + t + 1) L + u; Q+r (t 2 + u?t + u;) L + ^ Q+r \t 2 +ut + u 2 ) L ). 



1 "Superintegrable": if integrable both due to Onsager's algebra and a Yang-Baxter equation. 
2 For reviews see e.g. |l6|,[|17f 
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3 Z 2 : The Ising case 

In flU) we have seen three different variables in which to write the polynomials: z, c and s. 
Since for ( |i~3| ) we only need the zeros Zj or Cj or Sj, the choice which polynomials we should 
prefer can be decided by convenience. The range in which the zeros appear is dictated by the 
hermiticity of Tv- 8 ' which requires 

- 1 < Cj < +1 or - oo < Sj < or \zj\ = 1. (18) 

Let us warm up by considering the Ising case N = 2 for which there are two charge sectors 
Q = 0, 1 and we have s = t 2 . From fli~5| ) Baxter's polynomials are 

Y 



Pq\s = t 2 ) = {(t + 1) L + (-l) Q (t - 1) L } ; a L)Q = L + Q mod 2. (19) 



Their zeros are well known 11121 



•an 2 <! - ( j 



vr /. 1-Q 



One easily derives (see e.g. the next section) that for L even these polynomials satisfy the 
differential equation (DE) 

^s-l)^-2f Q ^ +9Q P^=0. 



with f = s(2L - 3) + 1; g = L(L - 1); 

h = s(2L - 5) + 3; 9l = (L - 1)(L - 2) (20) 

(for L odd interchange /q «-> /i and go ^ Si)- These are hypergeometric DEs with the polyno- 
mial solutions 

P^ L) ( S ) ~ 2 Fx(3±^, 8^; Q+\; s) for L even, 

4 L) ( S ) ~ 2^1 (i- 2 ^, i- 2 ^; i-<9; s ) for L odd - 

How does the same information appear in terms of the variable c (it is the Cj which are directly 
required in (p^))? We define polynomials I^q\c) by 

ITq L) ( C ) = (c+ l)fc.«pf > (a = ^) , (21) 

where bL,Q was given in (pi]). After some algebra, the DE ( |20| ) can be rewritten in terms of 
11^ (c) and we find 
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These are Chebyshev-type equations. Their polynomial solutions, correctly normalized in ac- 
cordance with (21), are: 



Ilf } (c) = 2 k T k (c); nf \c) = 2 k U k ^(c); n ( 2fc+1) (c) = 2 k (U k (c) ± U k ^(c)), 

i 

for k = 0, 1, 2, . . . and we use 

T n (x) = cos (n arccos x); U n (x) = sin{(n+ 1) arccos x}/\/l — x 2 . 

As classical orthogonal polynomials the Chebyshev polynomials satisfy three-term pure recursion 
relations. In terms of the IIq^ these become: 

n^ +4) -4 C n^ +2) + 4nS 3 i) =o. (23) 

A main concern later in this talk will be to investigate whether for N > 3 the polynomials TIq (c) 
still form orthogonal sequences or at least preserve some features of the Chebyshev polynomials. 
We also consider briefly eqs,(|l~9|) rewritten in terms of the variable z. The interpolation 



polynomials J 7 ^ (z) are related to the P^ (s) by 

^(z) = (z + l) 2b ^ P W = (|_i) 2 ) (24 ) 

and for N = 2 we get 

Leven: 2^ A L) = z L + 1; 2 1 ' 1 f[ L) = L odd: 2 1 " i ^ L) = 

z z — 1 iz=pl 

For L prime these are standard cyclotomic polynomials Q k (z), in some other cases products of 
Q k (z), e.g. T\ 2 {z) = 2 11 Q 3 Q A Q 6 Q l2 ; ^ 13 (z) = 2 12 Q 13 ; F\\z) = 2 12 Q 26 . The correspond- 
ing Z3 interpolating polynomials J-q'\z) are more complicated and for them we have not yet 
discovered particular interesting features. 



4 ^-Polynomials Pq (s) 

To derive recursion relations and differential equations it proves convenient to start with the 
Baxter polynomials P < q\s). We shall describe the details of the derivations for the case Z3. 
For higher N the same methods can be used, just the formulae become more involved. Only 
occasionally we will quote the result for general N. 



4.1 Recursion relations 

Our starting point is, that equivalent to Baxter's definition (|T5|), we can define the P%\s) via 
the expansion of the multinomial (t 1 + t N ~ 2 + . . . + t + 1) L and collecting the terms with 
powers of t a for each a mod iV = 0, 1, . . . , N — 1. So, for iV = 3we write 

{t 2 + t + l) L = t a ^P^ L \s)+t aL ^P ( l L) {s)+t aL ^P^ L \s), (25) 
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with cll,q = %(L + Q) mod 3 , and the P^ are required to depend on s = i 3 only. From (|2 
we get immediately 

(t 2 + t + = (t 2 + t + 1) (t°^°P (L) (s) + i^pf^s) + ^> 2 P 2 (L) (s)) 



(26) 



Collecting the coefficients of mod 3-powers of t, and choosing L such that L mod 3 = 0, this 
gives 



/ p (i+i) 



(L+l) 




( P {L) 
(L) 



( pf +2 > 



\p 



P 



(L+2) 





(L+2) 




P 



One rotates the column vectors once more to get the Pq + ^ in terms of the P). L+ ~' 



Q 



(L+l) 



(27) 



Still taking L mod 3 = 0, by repeated application of fl27|) etc. we obtain 4-term pure recursion 
relations between polynomials of the same Q, e.g.: 



p(L+3) 
M) 




-3sP (L+2) 


+ 3s(s 


-1) 


p(L- 
M) 


-1) 


+ (« 


"I) 2 


p(L) 
M) 


= 


p(L+3) 
r \ 




_ 3 p(L+2) 


+ 3(s 


- 1) 


p(L- 


U) 


+ (« 


"I) 2 


p(L) 


= 


p(L+3) 
r 2 




-3sP 2 (L+2) 


+ 3(s 


-1) 


p(LA 
r 2 


-1) 


+ (« 


"I) 2 


p(L) 
r 2 


= 


p(£+4) 
Q 


- (2s 


+l)p( L+3) 


+ 2(s 


-I) 2 


p(L- 
Q 


hi) 


+ (« 


-I) 3 


p(L) 
Q 


= 


„p(i+4) 
S-Tq 


- (2s 


+ 1) p (L+3) 


+ 2s(s 


-I) 2 


M) 


-1) 


+ (* 


"I) 3 


p(L) 



= 0. 



for Q = 1, 2 

(28) 

In each except the fourth equations pairs of polynomials appear which have the same degree. 
The following recursion is valid for all L > and Q: 



3 (L+9) 

Q 



3((s - l) 2 + 9s)P^ +6) + 3(s - l) 4 P^ L+3) - (s - l) 6 Pf } = 0. 



(29) 



4.2 Differential equations 

Differentiating (^), for L mod 3 = we get: 



(2 t + l)L( P (L) + t 2 P x (L) + tP 



(t 2 + t + 1) {3t 2 (P (i)/ + t 2 ^' + *P 2 (L) ') +2iP[ L) + P 2 (L) } 



where Pq"' = dPffl /ds, and comparing the coefficients of mod 3-powers of t: 



( P^' \ 



3s(s - 1) 



P] 
P 



(L)' 



-B 



( p{L) 

(L) 
L 

(L) 



P 

V p 



2sL — sL — 

13 = ( — L 2sL- 2(s - 1) — L 

-L -sL 2sL-(s-l) 



(30) 



Writing ( p0[) shorthand as 3s(s — 1)7-" = B V, and differentiating again we have 

9s 2 (s-l) 2 P" = {£ 2 -3(2s-l)S + 3s(s-l)£'}P 
27s 3 (s - 1) V" = {B 3 - 9(2s - 1)B 2 + 18 (3s(s - 1) + 1)B + 3s(s - 1)(2B' £ + BB') 

- 18 s(s -1) (2s -1)B'}V (31) 
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We can use these relations to obtain decoupled DEs for the three charge sectors as follows: We 
form the three expressions (Q = 0, 1, 2) 



27s 2 (s - l) 2 P { Q LY " - 27 s(s - 1) f Q P$>" + 3g Q P { Q L) ' - (L - 1) h Q P { Q L) = 0. 



(32) 



which contain 9 coefficients /q, gQ and hQ. Then using ( plj ) we express the P^' , P V q> and 

Pq^ in terms of Pq^ which gives a matrix equation A4V = 0. Requiring A4 to be diagonal 
gives 9 linear equations for the coefficients /q , <?q , /iq . Solving these, the result for L mod 3 = 
is 



s(2L - 4 - Q) + 2 + Q; 



Q = 3 — Q mod 3; 



/; 
/; 2 



3s(4s - 1)L 2 - 3sL |4(s - 1)Q + 10s - 7} 

+(s - 1) {3Q 2 (s - 1) + 3Q(5s - 1) + 20s - 2)} 



{(8s + l)L-4(s-l)}L; 



hi 



{(8s + 1)L — 12(s — 1)} (L — 2); 



3s + 1)L 2 - (16s - 7)L + 6(s - 1). 



(33) 



For the other values of L mod 3 we get the same expressions, just for rotated values of Q . E.g. 
we have the same equation for L mod 3 = 0, Q = as for L mod 3 = 1, Q = 2. These DE's 



(32) are not anti-selfadjoint as one would like them to be. 

For higher TV" the derivations are completely analogous, just quite lengthy. One obtains 
Nth-oider DE's for the Z^r-polynomials. For L mod N = the relation generalizing ( |30| ) is 

N a (a-l)¥r=BWv 
ds 



( M 


-sL 


-sL 


-sL 


-L 


N-(N-l)(s-l) 


-L 


-L 


-L 


-sL 


Af-(N-2)(s-l) 


-L 


-L 


-sL 


-sL 


M-(N-S)(s 


V ••• 









using the abbreviation Af = (N — l)s L = N shift with bi,Q being defined in ([16] 



(34) 



/ 



4.3 Approximation of the zeros 

In (2^] it was shown that changing the variable t into [3 according to 

sin (/?- jj) 



(35) 



sin (/?+§)' 

the zeros of the polynomials come approximately equidistant in the interval — j? < j3 < j? 
Without approximations, for Z3 the equation Vl,q{£) = becomes 



(2cos/3r i = 2cos(l{/? + ^} + ^Y 



(36) 



Neglecting the left-hand side of ( p6| ) for L ^ 1 (this is exponentially good for |/3| <C ^ and 



dangerous only for 2cos/3 « 1, i.e. « ±^ ) we get the approximate solutions 



R ~ kt 71 *r 6fc + 2Q - 3 
Pk ~ K --; K = vr; 



A; = 1,2, 



<bL,Q- 



(37) 
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Table 1: Z 3 -Polynomials ILg'(c) for L < 8. 



L 




n[ L) 






ni L) 




I 


1 


i 






I 




2 


3c -4- 1 


3c — 1 






Q 




q 
o 


yc — o 


ac + 6 






yc — o 




4 


27c + 18c — 5 


T7 2 

27c — 


11 




27c — 18c — 5 




■5 


81c 3 + 27c 2 -57c -11 


81c 3 - 


27c 2 - 


57c + 11 


81c 2 - 21 




6 


243c 4 - 270c 2 + 43 


243c 3 - 


f 81c 2 


- 135c- 21 


243c 3 - 81c 2 - 135c 


+ 21 


7 


3 6 c 4 +486c 3 - 549c 2 - 270c+43 


3 6 c 4 - 


702c 2 - 


f 85 


3 6 c 4 -486c 3 -549c 2 4 


-270c+43 


8 


3 7 c 5 + 3 6 c 4 - 2754c 3 - 702c 2 


3 7 c 5 - 


3 6 c 4 - 


2754c 3 + 702c 2 


3 7 c 4 - 1782c 2 + 171 






+711c+85 




+71 lc 


- 85 







or 

l + s k ^ sin 3 (K + f) -sin 3 K 
° k ~ 1-sjfe ~ ~sin 3 (iT + f) + sin 3 iT ^ ^ 

For N = 2 the analogous approach leads to exact results. 

Here we shall not further discuss this useful trigonometric mapping of the zeros to the interval 
— jr < (3 < jj. Rather we focus our attention to the rational mapping from — oo < s < to 
— 1 < c < 1 and study whether traces of the orthogonality present in the N = 2 case survive for 
N > 3. 



5 The polynomials Hq (c) 



Q 

"-Q 



The general definition of the IIq (c) has already been given in (21). For illustration Table 1 
lists the n^ ) for L < 8. The inversion relations ([l?]) simply imply that 

IL%\c) = (-l)^QU^_ Q _ L (-c), (39) 

from which it follows that for L + 2Q mod N = the polynomials IIq for N odd are functions 
of c 2 only. 

5.1 Recursion relations 

We first derive recursion relations for the case Z3. We rewrite the recursive relations obtained 
for the Pq (s), ( |27|) in terms of the polynomials IIq' (c), with the result (taking L mod 3 = 0): 

n (L+m) x ( n (L +l) 



n (L+i+l) = c (i) 



n^ +1) J \ n 2 



Hi 





(L+i) 



(L+i) 



9 



C (0) 



c_ 

c_ c_ 



C (l) 



c_ 
c_ 

1 



C (2) 



1 1 




(40) 



with c± = c± 1. Several interesting relations follow immediately from the particular structure 
of the matrices : Still taking L such that L mod 3 = 0, we have 



n 



(L+l) 





(L+2) 



n 



(L+l) 



2 

(L+2) 



2(nr +n 2 



(L), 



n, 





(L+3) 



r(^+l) 



r(^+3) 



21^ 



(L) 



(L) 
1 

n^-n^ = 2nj L+1) ; n^ L+3) -n^ L+3) = 2n^ L+2) , (4i) 



which shows that the 9 polynomials of the subset with L mod 3 = 0,1,2 and Q = 0, 1, 2 are not 
independent. There are also derivative relations between these polynomials, see later eqs.(48), 



The relations ([40] ) mix the three charge sectors. In order to find recursion relations which do 
not mix the charge sectors Q, we form the matrices C, C 2 , C 3 : 



C = C (2) C (D C (0) 



/n 
n 
Vn 



(3) 


(3) 
1 

(3) 



C_|_ fi- 
ll 



n 

(3)' 



(3) 



c _n 



(3) 



Cj-C_ 



•+n 2 

n (3) 



(3) 



n, 



(3) 



9c 2 



nf } = 3c + i 



/ 



n. 



(3) 



3c- 1. 



C 8 




C + C_/3 2 + C+C_/3 2 - 



C+/?2+ 
0: 2 



c+c_/3 3 + c+c_/3 3 _ 



c+/33^ 

"3 



OL2 - 
fo± 



3 5 c 

= 3 5 c 3 T 3 4 c 2 



4 -270c 2 + 43 
5„3^ q 4„2 _ 135c±2 l 



a 3 = 3 8 c 6 - 3 7 • 5c 4 + 3 4 • 59c 2 - 341 
Aj± = 9(3 6 c 5 T 3 5 c 4 - 3 4 • 10c 3 

±3 2 • 22c 2 + 177c T 19) 



The matrix equation C 3 + /C 2 + (?C + /iII = contains three independent equations and has 



the solution / = -27c 2 + 15 



-3ni 3) ; g 



relation (valid for all L > and all Q = 0, 1, 2 )0: 



48; h = —64, corresponding to the pure recursion 



n 



(L+9) 



(27c 2 



15<+ 6 ) 



48 n 



(L+3) 



64 n 



0. 



(42) 



The degrees of the polynomials appearing in this relation increase by two from the right to 
the left. If we consider Q = L mod 3, then we have only polynomials in c 2 , so that the 
degrees appearing are consecutive in powers of z = 9c 2 ("simple sets of polynomials") with 
integer coefficients. These recursion relations are of simplest type (like those for the Chebyshev 
polynomials), with the coefficients not depending on L. Compare also the same equation (|2T 



for the P. 



Q 



In order to form simple sets out of the Hq (c) with the same Q having the degree in c increasing 
in steps of one without repetitions, every third L must be left out. E.g. for Q = 0, we can 
form the sequences L = 0, 2, 3, 5, 6, 8, 9, . . . or L = 1, 2, 4, 5, 7, 8, ... , i.e. leaving out either all 
L mod 3 = 1 or all L mod 3 = 0. Analogous sequences can be formed for Q = 1 and for Q = 2. 

the relations ( [42] ) are valid also for 



If we define fl 
L = — 1, -2. 



(-1) 



n 



(-i) 



= 0; n 



(-2) 



1. 

4 ' 



n (-2) _ _1 



10 



Using again ( |40| ) and solving similar matrix equations, we find the following 4-term pure recursion 
relations: For all L > with L mod 3 = we get 



+ (1 



n (X+5+i) + (1 



0,1,2 with Q = A-j. 



Q 



3c) n 

for j 

9c)n^ +3+j) + 

for j = 0,1,2 with Q 



0; 

2 - 



j mod 3; 



i2c)i4 L+2+j) +sn( L+j) 



0; 



3 — j, 2 — j mod 3. 



(43) 



(44) 



ID- The 
L|D and 



Classical orthogonal polynomials must have 3-term recursion relations (see e.g. 
separation of the polynomials Pq L ^ for a fixed Q into two sets: the ones selected in 
the ones left out, reminds of Konhauser biorthogonal polynomials |24j] (KBOP). KBOPs have 
4-term pure recursion relations if the first set consists of the polynomials in c and second set of 
polynomials in c 2 . In some cases of ([43"|) this is what we just did, (the Hq' + ^ are functions of 



c 2 ), but not in all cases: for j = 2 in ( |4"3| ) H^q 1 are left out and these depend on c. 

We can apply the same technique to find recursion relations also for higher Z^v polynomials. 
Generally, we find N + 1-term pure recursion relations. For Z4 the result analogous to (E5) is 



n 



(L+16) 

Q 



4n (4) n (L+12) 



128(14c 2 



17)lFg +8) 



2048cII^ +4) + 4096IlS 3 L) 



0. 



(45) 



with n 



(i) 




64c 3 - 56c. For Z fi 



n 



(L+25) 



5n< 5) ng- +20) + a 2 U^> + a 3 n^"> + a 4 n^ + a 5 n 



AL+15) 



(£+10) 



0. 



(46) 



with -5n[, 5) = -3125 c 4 + 3750 c 2 - 705; a 2 = -2 6 5 (375 c 2 - 383); 

a 3 = -2 10 (5 4 c 2 -45- 13); a 4 = 2 16 • 5; a 5 = -2 20 . 



There are also N + 1-term recursion relations analogous to (f43|), ( |44| ) for sequences formed of 
H^q\c) with consecutive degrees in c. Generally, for this, certain every Nth L must be left 
out. If these were KBOPs, for having N + 1-term pure recursions, there should be subsets of 
polynomials depending on c . However, this is not the case: for any N there are only subsets 
depending on c 2 because of (p9|). 

One might think that orthogonal polynomials which usefully approximate the could be 



obtained recursively, neglecting the last term in (42). Using monic polynomials fig, this looks 
promising because eq.fl42]) becomes 

g(L + 9) = ( c 2 _ 5\ fi (L+6) _ ^L + 3) _^_ fi W 

Q V 9/ Q 243 Q 19683 Q 

with the last coefficient looking small. However, new polynomials defined through (consid- 
ering L mod 3 = 0, Q = 0): 

Ri L+&) = ( c 2 - f| Rf +3) - ^R { . L) 



\ 9J 243" 

with = 3; R^ = c 2 — | don't seem to be useful, since the R^ have zeros which for 

increasing L spread out of the interval — 1 < c < +1, violating the hermiticity of the hamiltonian. 
The method using eq.(p6|) does not suffer from this deficiency. 
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5.2 Differential equations for the Uq\c) 



The derivative relation (Q) can be rewritten for the IIq , giving e.g. for L mod N = 0: 

/ no' \ 

c+IV 

c+n 2 ' 

c+IV 



V 



1 



/ 








Lc_ 


Lc- 




" ^ 


/ 


n 






Lc + 


-Nc + + 2 


Lc + 


Lc + 






c+Iii 






Lc + 


Lc- 


-Nc++A 


Lc + 






c+n 2 






Lc + 


Lc_ 


Lc- 


-iVc++6 






c+n 3 




V 












V 




) 



(47) 



where the primes denote differentiation d/dc and we omitted the superscripts (L) on the IIq. 
The first component of this equation gives the nice simple relation 

AT-l 



IV 



L 

N 



£n Q . 



(48) 



We now concentrate on the case iV = 3. From relations analogous to ( |47| ) we get the following 
simple derivative relations: 

3c_n ' = L(n x + n 2 ) - 2n ; 3c + n 2 ' = L(n + n x ) - 2n 2 

3(c 2 - i)ni' = L(c+n 2 + c _n ) - 2cUi. (49) 

Using these and similar relations, after considerable algebra, we get the DE 



27(c 2 -l) { (c 2 -l) 



dc 3 



a(c)- 



d^ 



(L) 



dc 2 



3{3(c 2 -l)L(L + l) + /3(c)}^_ 



(L) 



0. 



(50) 



We collect the coefficients in Table 2. Observe that these DE for the IIq (c) are much simpler 
Table 2: Coefficients appearing in the DE (pOt) for the Z 3 -polynomials U^\c) 



L mod 3 


Q 


a(c) 


/3(c) 


l(L,c) 








4c 


-18c 2 + 10 


-2cL 2 (L + 3) 





i 

2 


7c ± 1 


-90c 2 T 24c + 34 


-2cL 2 (L + 3) - 3(3c ± 1){L{L + 1) - 6) 


1 


1 


6c 


-60c 2 + 28 


-2c{L 2 (L + 6) + 3L-10} 


1 


2 



6c + 2 


-60c 2 ± 36c + 16 


-2c {L 2 (L + 6) + 3L - 10} ± 6(L(L + 1) -2) 


2 


2 


8c 


-126c 2 + 46 


-2c{L 2 (L + 9) + 6L-56} 


2 




l 


5c ± 1 


-36c 2 T 12c + 16 


-c {L 2 (2L + 9) + 3L - 4} T 3L(L + 1) 



than the corresponding eqs.(|32"|) for the Pq (s): there the coefficients /q of the second derivatives 

Pq 1 ^ p^ ) were L-dependent, which here is not the case. However, although in ( |50| ) the third 
and second derivative terms can be combined into a total derivative, we don't see how to put 
the whole DE into an anti-self conjugate form. 
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5.3 Separation property of the zeros. 

We want to collect further evidence that the polynomials IIq violate features which are crucial 

for classical orthogonal polynomials. Table 3 shows the zeros of the polynomials H-J^ for low L. 
Fully satisfied is the confinement of the zeros Cj to the interior of the basic interval — 1 < c,- < 
+1 (this must be so because of the hermiticity of the hamiltonian) . However, the separation 
property of the zeros of polynomials sequences successive in c can not be satisfied: If we take 



out the L = 3, 6, 9, . . .-polynomials as we did in (43) for j = 1 with Q = 4 — j then the zeros 



.78444; .77622 for L = 5,7 are in wrong order. Alternatively, taking out L = 4,7,11,... as 



in (43) for j = 2 with Q = 2 — j then .95835; .93330 for L = 6,8 are in wrong order too. 



Considering the polynomials in c 2 only (Q = and L = 3, 6, 9, . . . as in (|42|)) gives the correct 

,2 

r 



separation of the c 2 



Table 3: Zeros Cj of the Z3-polynoms ITq \c) showing the separation property violation at the 
upper corner near c = 1. 

_± Cj (j = 1, [2L/3]) 

2 -.33333 

3 -.74536 .74536 

4 -.87766 .21100 

5 -.93203 -.18575 .78444 

6 -.95835 -.43894 .43894 .95835 

7 -.97264 -.60038 .13013 .77622 

8 -.98106 -.70632 -.11267 .53341 .93330 

9 -.98634 -.77818 -.29702 .29702 .77818 .98634 

10 -.98983 -.82848 -.43636 .09047 .58910 .90843 

11 -.99222 -.86467 -.54241 -.08238 .39948 .77767 .97120 

12 -.99392 -.89136 -.62401 -.22459 .22459 .62401 .89136 .99392 



5.4 The Uq expressed in terms of Jacobi polynomials: 
Partial orthogonality 

In order to find out whether the (c) satisfy some generalized kind of orthogonality, we looked 
for the eventual weight functions. The definition interval is certainly — 1 < c < +1 and as a 
simple guess one may try Jacobi weights 

w«m = (i - c ) a (i+cy 3 . 

Since for Z2 the IIq'' 1 are Chebyshev polynomials, i.e. there a = f3 = ±l/iV, we tried whether 
e.g. the Z3 polynomials are related to a, [3 = ±| ± |-Jacobi polynomials. Indeed, expanding 

the Uq (c) in terms of Jacobi polynomials Pn°'^ (c) and playing around with several choices of 
a and b in a = a/N and f3 = b/N we find remarkable features which we report now. 
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Table 4: Expansion vectors ( ]5l| ) of the Z3-polynomials 6 \ 2L I^ IIq^ in terms of Jacobi polyno- 
mials. Note the many vanishing components ao, a%, etc. 



6 -[2L/3] n W( c ) for Z 3 


L 


Q = 


Q = 1 


Q = 2 



2 

3 

4 
5 
6 
7 
8 
9 
10 



0.3 



l l 



2 _1 
" 3 ' 3 



1 _1 
3' 3 



0' 9' 3 

[0,0, |] 2 _2 

3' 3 

[n _J_ o — 

L u ' 18 ' u ' 20 
L u ' u ' 45' 35' 35 
L u ' u ' 36 ' u ' 70 J 

[0,0,0,-^,0,^], 
[0 ^ ^ ^ 

L u ' u ' u ' 756' 105' 154' 

[0,0,0,0,-^,0, 



2 _1 
3' 3 



2 _2 

3 ' 3 



0,| 



1 1 

"3'3 



1 1 
"3'3 



1 1 
"3'3 



[°'i]i_i 

3' 3 
L u ' 4' 2 
L u ' 18' u ' 20 

[0,0,0, |> ]i _ a 

u 3 ' 3 

[o ^ ^ I 

L u ' u ' 63 ' 40' 70 J 

[0,0,0,-^,0,^ 

[0, 0, 0, 3go, 0, 

[n — — - 

|_u,u, u, 360 , 11Q , 56 , 30g 



i i 

3 ■ 3 



1 1 
3'3 



I 2 

3 ■ 3 



1 1 
"3'3 



2 2 
"3'3 



1 2 

3'3 



[o,|L 

[0,0,i]_ 

\0 i -5-1 

L u ' 30' 20 J 

[0,0,0,^ 
[0,0 
[0,0 
[o, 0, 0, ggp, 0, 
[0,0,0,0,^,0, 



1 1 

3'3 



36 ' u ' 70 



111 

168 ' 105 ' 28 



1 2 

3'3 



1 1 

"3'3 



9 
308 



2 2 
"3'3 



We use the standard definitions given e.g. in Szego |2l| with 

P (^) (1) = + P («-«( C ) = 1; PfHc) = (-IT P^ a H-c). 

We present our results in Table 4 using the notation 

n 

E^P^c). (51) 



ao, ai, . . . , a n 



k=0 



(n can be seen from the number of arguments). Surprisingly, when making a convenient chioce 
of a and /?, with increasing L an increasing number of the lowest components of the vectors 
[ao, ai, . . turns out to be exactly zero. So, indeed our ^-polynomials are related to Jacobi 
p( a >@\ Defining a scalar product 



(n£\ n%X P = J 1 dc (i - c r(i + c)/ 3 njf (c) n$\c) 

by inspection of Table 4 for JV = 3, we conjecture the "partial orthogonality" 



(nf >, n 3fc,) ) 



2_i=0 for 2fc + 1 < fc'. 

3 ' 3 



(52) 



(53) 



because the expansion vector of IIq^ has 2k + 1 components and for k < 35 we have checked that 
its k lowest components are zero. Interestingly, the first non-zero components are very small 
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Table 5: Expansion (|54j) of Z4-polynomials 2 [ 3i / 4 l Hq in terms of Chebyshev polynomials. 



The same components appear (partly shifted) e.g. for 7q+i, 72+3 and 81+3. 



T 

L Qi+Q 2 


0-[3i/4] 1 m( L )f^i _L rr^)^^ 


J 


0-[3i/4] TT^V^ 
Z 1 ' J ILq [Cj 


ll+2 


r-i i 


2i 


[0,2] T 


^0+2 


Tn 1 1 


23 


[ 2 \u 


3o+i 


[0,0, 2] T 


4 


[0,-l,0,2] T 


Q 


[°> 2 Jt/ 


A 

42 


r 1 n ol 


A 

4l+3 


[0,0,% 


5o+3 


Tn 1 r\ /il 

[0, 1,0,4]^ 


5l+2 


[0,-3,0,4]^ 


61 


[0,0,-2,0,8] T 


DO+2 


[0,0, -3,0,4]^ 


03 


\c\ n ol 

[0,-2,0,% 


7o+i 


[0,0,0,-4,0,8] T 


80 


[0,0,1,0,-8,0,8] T 


'2+3 


[0,0, -4,0,8]^ 


82 


Tn 1 n q n Ql 
[0, 1,0, -8,0,% 


O 

oi+3 


[0,0,0, -4,0, S\ u 


»0+3 


[n n on A Ci 1 £il 

[0,0, -3,0, -4,0, lbjy 


9l+2 


[0,0,5,0, -20,0, 16] j; 


10i 


[0,0,0, 2,0, —24,0, 32] T 


10o+2 


[0,0,0,5,0,-20,0, 16] v 


10 3 


[0,0,2,0,-24,0,3% 


H0+1 


[0, 0,0, 0,6, 0,-32,0, 32] T 


12o 


[0,0,0,-1,0, 18, 0,-48, 0,32] T 


H2+3 


[0, 0, 0,6, 0,-32,0, 32] u 


12 2 


[0, 0,-1, 0,18, 0,-48,0, 32] „ 


12l+3 


[0, 0, 0,0, 6, 0,-32,0, 32] ^ 


13o+3 


[0,0,0,5,0,-8,0,-48,0,6% 



and increase strongly before reaching a plateau level around the kth component: E.g. for IZq 
we have 00 = . . . = 034 = 0; 035 = —0.0007395; 0,3$ = 14.6665; 052 = —3.186 • 10 9 . It is a simple 
consequence of the reflection relation (|3^ ) that for fixed L always one pair of Z3 charge sectors 
differs only by the interchange of a and (5. 

Several other partial orthogonality relations can be inferred from other vanishing first compo- 
nents in Table 4. 

For the ^-polynomials we find that these show a similar partial orthogonality if expanded in 
terms of Chebyshev polynomials U n (c) and T n (c), i.e. for a = (3 = ±i. Here using the normali- 

(±1+1) 

sation of the Chebyshev polynomials leads to much simpler coefficients than using P n 2 ' 2 (c) , 
so we define 

n 

[op, qi, ...,an)u =^ dkUkjc), (54) 

fc=0 

analogously for T. Table 5 lists results for L < 13. Observe the remarkable symmetry between 
pairs of sectors obtained by replacing Tp. by U^-i (e.g. for n 1 ^ 1 ^ and IIq 1 ^ 2 ^ taking 
L mod 4 = 0) or, replacing Tfc by Uk for and n 1 ^ 1 ^ (taking still L mod 4 = 0). These 

relations are not at all simple if expressed in terms of the 11^ (c). 
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We conclude these observations giving in Table 6 few results for N = 5 and N = 6: For 
N = 5 we find partial orthogonality due to vanishing first components with a = —j3 = =t| and 
a = —(3 = ±g according to the sectors L mod 5 and Q considered. For N = 6 both Chebyshev- 
and e.g. a = —(3 = i^-Jacobi weights are useful. As for N = 4, various symmetries show up, 
e.g. for N = 6, L mod 6 = exchanging Q = with Q = 3 corresponds to exchanging the 
Chebyshev polynomials and Uk-i- A convenient method to perform expansions in terms of 
Jacobi polynomials is given in |p5| , Sec. 7.2. 



6 Conclusions 

We have studied various properties of the polynomials, which via their zeros determine the energy 
eigenvalues of hamiltonians satisfying Onsager's algebra. Three different versions have been 
considered with the zeros in the hermitian case on the negative real axis, the interval (—1, +1), 
and the unit circle, respectively. The most simple properties emerged in the second case, the 
polynomials Hq'\c) which for Z2 are Chebyshev polynomials. For the general Zjv we argue 
that the polynomials considered have no shorter than N + 1-term pure recursion relations. This 
excludes that the Ilg are classical orthogonal polynomials. Konhauser biorthogonality (which 

allows 4-term and higher recursion relations) is not seen to be a property of the Hq . However, 
the expansion in terms of Jacobi polynomials reveals a very remarkable partial orthogonality with 
respect to Jacobi weight functions. The deeper meaning of the latter remains to be studied. 



There are more general definitions of biorthogonal polynomials in the literature [26], e.g. of 



Iserles and N0rsett p7l] and Van Iseghem [28]. For hypergeometric functions 2F2 there are four- 



term recursion relations, see e.g. [|29[] Chap. 14. However, these have a much more complicated 
structure than the recursions which we foun 
possible relevance of these latter structures. 



structure than the recursions which we found for the More work is needed to clarify the 
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Table 6: Expansion vectors for some Z5- and Zg-polynomials. Since in several cases we don't 
show all components but just the first and last ones, in the column c£zm[] we give the lengths of 
the component vectors. 



L Q dim[] 5"[ L / 5 1 Ill L) (c) for Z 5 



4 2 
5i 
61 

60 

3 

9 2 
14 2 

192 

24 2 



3 
4 
5 
5 

7 
11 

15 

19 



[o,o,fM f ,- f 

' 

[0,-8,0,5013 _s 

5 ' 5 

[o,o,=fe,o,^] f) 

[0,^16,0, ±150, ^ 



5 • 5 



Tn n 4400 n - 2 60000 n i 2 50000 



0,0,0,0 



77 

-179 2 Q00 
143 

- 2 Q48000 



231 



0^0 



2 _ 2 

5 ' 5 

9 2 -10 8 







1 2 5-10 8 



,0, 



-,o,...,o, 



0,0,0,0, 221 > u > 2 1 
0,0,0,0,0,0, ^»,0,...,0, 



561 24453 ' u ' 46189 

128-10 5 



-55-10 12 n 25-10 11 
2028117 ' u > 200583 



-688- 10 15 
420756273 



,o, 



5-10 16 
90751353 



3 3 
8 3 

13 3 
I83 
38 3 



2 
6 
10 

14 
30 



[0,25 

[0,0,0, 



5 ' 5 

-2000 



0. 



125000 
63 



10 8 q 25-10 7 



000 - 9600 74 ' 1QJ 

' u ' u ' 11 ' u ' 273 ' ' 7293 ' ' 2431 

00000 1056 105 25 - 1Ql ° 



5 2 003 



1 _i 

5 ' 5 



nnnnnnnnn 478150656 - 10 o n 625,10 5 

u, u, u, u, u, u, u, u, u, 2 68801 ' u ' " " " ' u > 375840831 2 44 2 63 





dim[ ] 


6 -[L/6] n (i)( c j for Ze 




3 4 


2 


[0,36] 11 

6 1 6 




4 3 


3 


[0, 0, 1441 1 1 

3' 3 




10 3 


8 


[0,0,0,8736,0, - 50 9 4 f 60 ,0, ^Hf 2 ], 1 




60 


6 


[0,f ,0,-90, 0,81] T 




63 


5 


[f ,0,-90, 0,81] ^ 
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